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Abstract 

If interpreted in terms of decaying dark matter, the steep rise in 
the positron fraction of cosmic rays above 10 GeV, as observed by the 
PAMELA experiment, suggests an underlying production mechanism 
that favors leptonic channels. We consider a scenario where a portion of 
the dark matter is made of the gauginos of an unbroken hidden-sector 
U{l)x, which interact with the visible sector only through a tiny ki- 
netic mixing. The second component of the dark matter is made of 
neutralinos, and depending on the mass spectrum, the lightest neu- 
tralino or the hidden gaugino becomes unstable and subject to decay. 
We analyze the cosmic rays, namely the contributions to the positron, 
the extragalactic gamma-ray and the antiproton flux, which poten- 
tially result from these decays and demonstrate that the production 
of antiprotons can be naturally suppressed. Furthermore, we briefly 
discuss the apparent double-peak structure of the ATIC data in light 
of cascade-decaying hidden gauginos, as well as possible signatures at 
Fermi. 



1 Introduction 



The PAMELA collaboration has recently reported a measurement of the 
positron fraction [I] which shows a significant excess at high energies com- 
pared to the expectations from spallation of cosmic rays on the interstellar 
medium [2]. This result confirms the excess in the positron fraction reported 
in the past by several experiments: HEAT 94/95/00 [3j, CAPRICE94 [3] and 
AMS-01 More importantly, the PAMELA collaboration has provided 
an accurate measurement of the energy spectrum of the positron fraction, 
revealing a steep rise between 7 and 100 GeV, possibly extending to even 
higher energies. 

The excess in the positron fraction has been interpreted as an indication 
for the annihilation [U [71 [U [9] or the decay [TTl [12] of dark matter par- 
ticles. If interpreted as dark matter annihilations, the dark matter particle 
must annihilate into W, Z, e, //, r for any dark matter mass, or into q, b, 
t, h for a multi-TeV dark matter mass. In both cases one needs large boost 
factors, e.g. from dark matter substructures, since the required annihila- 
tion rates are much larger than the cross-section suggested by cosmology, 
(av) ^3x 10-26 cmVs [131, would lead one to expect. On the other hand, 
if interpreted as dark matter decay, the dark matter particles must have a 
mass larger than ~ 300 GeV, a lifetime around lO^^s and the decay must 
proceed preferentially into charged leptons of the first and second genera- 
tion [l^. The properties of the dark matter particles are also subject to 
constraints from the diffuse gamma-ray flux as measured by the EGRET 
instrument [H], and the antiproton flux, as measured by PAMELA |16j . 
BESS95 dZ], BESS95/97 [E], CAPRICE94 [19j, CAPRICE98 [20] and 
IMAX [21]. More concretely, the good agreement of the theoretical pre- 
dictions for the antiproton flux with the measurements suggests that the 
dark matter particle annihilates or decays mostly into leptons. Moreover, 
taking additionally into account the excess in the total electron+positron 
flux as observed by the ATIC experiment [22] hints to dark matter particle 
masses above 1 TeV0 Although the interpretation of the positron excess in 
terms of dark matter is very suggestive, it should be borne in mind that 
nearby astrophysical sources such as pulsars might produce sizable positron 
fluxes in the energy range explored by PAMELA [28j. 

In this paper we will discuss the cosmic-ray signatures of super symmetric 
scenarios with a gaugino of an unbroken hidden-sector U{1) gauge group 

^For the combined electron+positron flux, we use results from measurements that have 
been undertaken by ATIC, HEAT94/95 [23], PPB-BETS [24], H.E.S.S. [25], BETS [26], 
CAPRICE94 4 and SANRIKU [27]. 
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which interacts with the visible sector only via a tiny kinetic mixing [29] . 
We will speculate that the tiny kinetic mixing induces the decay of the 
dark matter particles into lighter supersymmetric particles and positrons, 
thus providing a potential explanation to the excess observed by PAMELA. 
Assuming exact -R-parity conservation, two possibilities may arise. First, we 
will study the case that the hidden gaugino mass is smaller than the lightest 
neutralino mass. We will show that the neutralino mostly decays into two 
charged leptons and the hidden gaugino. The hard positrons produced in 
the decay can then potentially explain the steep rise in the positron fraction 
observed by PAMELA and the absence of an excess in the antiproton flux. 
Secondly, we will study the case that the lightest neutralino mass is smaller 
than the hidden gaugino mass. If this is the case, it is the hidden gaugino 
which decays into the lightest neutralino, either directly or in a cascade 
decay, when there are supersymmetric particles with masses between the 
hidden gaugino mass and the lightest neutralino mass. Either particle could 
be the dominant component of dark matter. However, we will focus in this 
work on the possibility that the dominant component of dark matter is the 
lightest neutralino, which may allow direct dark matter detection. 

In section 2 we will shortly review the main features and possible origins 
of scenarios with a hidden abelian gauge group and kinetic mixing. In 
section 3 we will review the procedure to calculate the gamma-ray, positron 
and antiproton fluxes at Earth from dark matter decay. In section 4 we will 
present our results for the decaying neutralino case and the decaying hidden 
gaugino case. Lastly, in section 5 we will present our conclusions. 



Many extensions of the Minimal Supersymmetric Standard Model (MSSM) 
contemplate the possibility of a hidden sector, consisting of superfields which 
are singlets under the Standard Model gauge group. Hidden sector super- 
fields usually couple very weakly to our observable sector, thus constituting 
a very natural arena for finding dark matter candidates. We consider an 
extension of the MSSM by a hidden abelian gauge group U{l)x (for details 
about the model see Ref. [29]). This gauge group remains unbroken at low 
energies and couples to the MSSM only through a tiny kinetic mixing, x, 
with the hypercharge U{1)y |30j : 



2 Model 





kin. mixing x 
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We assume that all matter states charged under U{l)x are heavy and negli- 
gible. As a consequence the U{l)x gauge boson completely decouples from 
the observable sector. However, as shown in Ref. [29], a non-zero mass- 
mixing of the order of 6M ~ 0(x • Mx) between hidden gaugino and bino 
generally remains. More precisely, in the basis where the kinetic terms are 
canonical, the extended (5 x 5) neutralino mass matrix reads, to lowest order 
in X, 



(Mx 


6M 








\ 


6M 


Ml 





-MzCj^sw 


Mzspsw 








M2 


Mzcpcw 


-Mzspcw 





-MzC/3Sw 


Mzcpcw 







V 


Mzspsw 


-Mzspcw 




/ 



Here, /i denotes the MSSM |U-term, Mz the mass of the Z° gauge boson. Mi 
and M2 the bino and wino masses, respectively, sw the sine of the Weinberg 
angle and sp \s related to the ratio of the two Higgs VEVs. 

Below, we will concentrate on the case where the lightest super symmetric 
particle in the visible sector is a bino-like neutralino Xi- Then, depending 
on the masses of the hidden gaugino and of the neutralino, one of the two 
particles becomes unstable with a lifetime that is roughly given by 

..,,~O(10--10)xl„-..(i|^)"(^)", (3) 

where we made use of the mixing angle 9 5M/\M^() — Mx\ ~ C'(x)- The 
exact prefactor depends on the dominant decay modes and on the mass 
spectrum of the supersymmetric particles. However, it is apparent that a 
lifetime around 10^^ s, as required to fit the PAMELA excess with decaying 
dark matter particles, implies an extremely small mixing of the order of x ~ 
10-24. Note that the mixing must be somewhat larger when the decaying 
particle is only a subdominant component of the dark matter, but an upper 
bound of roughly ~ 10^^'^ holds from the requirement that the particle has 
not already decayed. 

The thermal production of hidden gauginos by oscillations between bino 
and hidden gaugino, which generally takes place in the primeval MSSM 
plasma, is irrelevant for the mixing parameters that we are looking at ^29j . 
However, the hidden gaugino may be produced non-thermally, e.g. in the 
decay of a heavy gravitino, or it may be a thermal relic of the hidden sec- 
tor [31] . In the latter case one requires additional particles that are charged 
under the hidden U{l)x with masses around Mx, however. For simplicity. 
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we will assume throughout this work that the lightest neutralino in the vis- 
ible sector has the right energy density to make up the dominant part of 
the observed dark matter, whereas the abundance of the hidden gaugino is 
always sub dominant, px <C p^o ~ pbm- 

The actual mass scale of the hidden gaugino depends on how the breaking 
of supersymmetry is mediated to the visible sector and the hidden U{l)x- 
If the soft masses ~ Msoft in the visible sector arise from gauge mediation, 
whereas U{l)x couples to the supersymmetry-breaking sector only gravita- 
tionally, the predicted mass hierarchy is Mgoft ^ ^-3/2 ~ where 
denotes the mass of the gravitino. In this scenario the lightest neutralino 
cannot be dark matter, because it would decay into the gravitino in the early 
Universe. However, if we suppose gravity mediation to the visible sector and 
anomaly mediation [32] to the U{l)x, one expects Msoft ~ ^-3/2 ^ Mx- If 
the soft masses of both the visible sector and the U{l)x arise from grav- 
ity mediation, one in general expects that all masses are of the same order 
-Msoft ~ ^3/2 ~ Mx- Below we will assume that the gravitino is heavy 
enough to have no impact on the decay modes, in^/2 > max(Mx, M^o). 

Additional hidden sector U{1) gauge factors are a generic feature of 
string compactifications. For example, in the "mini-landscape" of orbifold 
compactifications of the heterotic string [33] one encounters, at the com- 
pactification scale, a breaking of the gauge symmetry to a theory involving 
many hidden ;7(l)s, e.g. Eg x Eg Gsm x U(1)^ x [S0(8) x SU(2) x U(1)^] 
and the like. Similarly, type II compactifications generically invoke hidden 
sector [/(l)s, often also for global consistency requirements. Some of these 
hidden U{l)s may remain unbroken down to very small scales |i34j. 

Kinetic mixing is generated by the exchange of heavy messengers that 
couple both to the hypercharge U{1)y as well as to the hidden U{l)x- 
Correspondingly, it is loop suppressed, x = dYdx C/{167r'^), where gy and 
gx are the abelian gauge couplings and C is a dimensionless constant. In 
field-theoretic setups, the latter is naturally of order one [30] and thus way 
too large for our purposes. However, it can be much smaller if there are 
additional gauge or global symmetries (cf. Refs. [29^ I35j). Moreover, in 
models arising from string compactifications rather small mixings seem to 
be generic [36l[37l[34l[38]. In the context of compactifications of the heterotic 
string, the size of the mixing has been estimated as |36j 

gygx Am 
^ " 167r2 Mp ' 

where Am <C Mp is the mass splitting in the messenger sector and Mp 



5 



the Planck scale. Indeed, the mixing is quite small in this case {e.g. x ~ 
10^^^^, if Am is associated with gauge mediated supersymmetry breaking, 
Am ~ 100 TeV), but still to large for our purposes. However, a sufficiently 
strong suppression of the coupling could be achieved in models with multiple 
hidden U{l)s t29j. 

In compactifications of type-II string theories, where the hidden ?7(l)s 
arise via D-branes in the (extra-dimensional) bulk that have no intersection 
with the branes responsible for the visible sector, one has to distinguish two 
cases. 

i) In large volume scenarios, where the hidden D3-branes and the branes 
supporting hypercharge are located at generic positions in some Calabi-Yau, 
kinetic mixing is suppressed by the large volume Vq of the extra dimen- 
sions [38] , 

X.^If^(ViMf)-", (5) 

where Mg is the string scale, which is related to the Planck scale via Mg ~ 
QsMp/ y^VgMp, where Qs is the string coupling. We thus see that we may 
obtain a kinetic mixing of order 10~^^, if we assume V^M^ ~ 10^^, corre- 
sponding to a string scale of order Mg ~ TeV. However, this case is prob- 
lematic because the mass of the gravitino is of the order of meV [39j, with 
impact on the stability of super symmetric dark matter candidates like the 
neutralino or hidden gaugino. 

ii) In scenarios with significant warping, such as KKLT |40j, the standard 
model stack of branes, notably the brane featuring the hypercharge U{1), 
is placed at a special position - at the tip of a warped throat - while the 
hidden brane is separated from it by a distance d along the throat. In this 
case kinetic mixing may be exponentially suppressed [38] like 

(6) 

reminiscent of a "Yukawa type" interaction. In fact, in warped compactifi- 
cations, the closed string fields that mediate the kinetic mixing, notably the 
NS-NS B form and the R-R C form fields, acquire masses, 

m^nMs/{V^Ml), (7) 

for some integer n, from the vacuum expectation values of their three-form 
field strengths, H^^ = di?2, -^3 = dC2, which are threading three-cycles of 
volume V3 and cause the warping. Therefore, we conclude that no strong 
fine-tuning is required to obtain a tiny kinetic mixing in a KKLT-like sce- 
nario. Indeed, the estimate ([6]) yields the desired small value x ~ 10~^^ 
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for md ~ 48. 



As mentioned above, depending on the mass spectrum of the supersym- 
metric particles, either the hghtest neutrahno, Xi; or the hidden gaugino, X, 
becomes unstable. The relevant decay modes are shown in Tab.[TJ These de- 
cays may be detectable as anomalous contributions to the cosmic-ray fluxes 
observed at Earth. We will discuss this in some detail in the next section. 



Neutralino / Hidden gaugino decay modes 



M^o > Mx 



X? 



ffl/R - ffX 

Xh' 

xz^ 



X 



M^o < Mx 

f f _ 

■I J L/R 

x°^° 
I xfw^ 



f fXi 



Table 1: Dominant decay modes. Depending on the masses of the hidden 
gaugino, Mx, and the lightest neutralino, M^o, one of the two particles 
becomes unstable with a lifetime roughly given by Eq. ([3D. Since the three- 
body decay into fermion pairs // is mainly mediated by virtual sfermions, 
/*, we show this explicitly. Furthermore, when a sfermion is lighter than 
the decaying particle, the corresponding three-body decay crosses over to a 
cascade decay. The subsequent decay and fragmentation of the Higgs- and 
gauge bosons, charginos and neutralinos is not shown. Note that the letter 
/ represents any lepton or quark. 



3 Propagation Models 

The decay of dark matter particles induced by the kinetic mixing discussed 
above will generate potentially sizable numbers of Standard Model particles 
which may be detectable in cosmic-ray experiments if the dark matter decay 
rate is high enough. In the following, we will consider the cosmic-ray signa- 
tures resulting from the decay of neutralinos and hidden gauginos in three 
different channels, namely gamma rays, positrons and antiproton^. Inter- 

■^Dark matter decay also produces a flux of neutrinos which are unfortunately unob- 
servable due to the large atmospheric neutrino backgrounds [41j . 
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estingly, in all of these channels, experiments are now reaching sensitivities 
to dark matter lifetimes of the order of 10^^ s. 



Gamma rays. Photons from the decay of GeV - TeV mass dark mat- 
ter particles will manifest themselves as an anomalous contribution to the 
diffuse extragalactic gamma-ray background. More specifically, there are 
two distinct contributions to the diffuse gamma-ray flux from dark matter 
decay. There is an anisotropic component from the decay of dark matter 
particles in the Milky Way halo, as well as an isotropic diffuse emission from 
dark matter decaying at cosmological distances. The latter contribution is 
red-shifted due to the expansion of the Universe. Thus, for a dark matter 
mass TTT-DM and a dark matter lifetime of tdmj the gamma-ray flux will be 
given by the sum 
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dE 



DM 



'2 dJ 
^ dE 



+ 



halo 



'2 dJ 
^ dE 



where the halo component is given by 



'2 dJ 
^ dE 



halo 



2^2 I r 

mDM dE Stttdm Jios 



PDM{l)dl , 



(8) 



(9) 



in which poM is the dark matter halo profile, which in the following we 
assume to be an Navarro- Frenk- White (NFW) [l3] profile for definiteness. 
The results will have only a slight dependence on the particular choice of 
halo profile. dN^/dE denotes the energy spectrum of photons from dark 
matter decay. The integration extends along the line of sight (los), yielding 
an angular dependence of the signal. For our results, we will average the 
anisotropic halo contribution over the whole sky excluding a band of ±10° 
around the Galactic disk. The extragalactic component is given by 



E' 



.dJ 
dE 



extra 



with the coefficient 



2E^ 



Cry 



dy 



dN^ 



,-3/2 



d{yE) yrriw^W^ 



(10) 



(11) 



The integration is over the red-shift z, with y = z + 1. In these expressions, 
r^A — SOm, — 0.24 and JIdm — 0.20 are the vacuum, matter and dark 
matter density parameters, respectively, while Hq = 70.1 kms~^ Mpc~^ is 
the present value of the Hubble parameter and pc is the critical density. One 
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finds that the halo component is dominant, although both contributions are 
of the same order. For the background to the dark matter signal from 
conventional astrophysical sources, we will assume a power law spectrum 
that fits the data points below 1 GeV, 

= 6.8 X IQ-'^E-^-^^ GeV (cm^ str s)"^ . (12) 

background 

Antimatter propagation. At the energies of interest here, only antipar- 
ticles created within the Milky Way's dark matter halo will be of importance. 
Antimatter propagation in the Galaxy is a fairly complicated process due 
to diffusion, energy loss and annihilation effects. It is commonly described 
using a stationary two-zone diffusion model with cylindrical boundary con- 
ditions [H]. Under this approximation, the number density of antiparticles 
per unit kinetic energy, /(T, r, t), satisfies the following transport equation, 
which applies both for positrons and antiprotons: 

^ = V-[K{T,r)Vf]+ (13) 
+^[6(r, f)f] - V • [V,{r}f] - 2h5{z)T^^^f + Q(r, r) . 

We assume free escape conditions for the solution /(T, f, t) at the boundary 
of the diffusion zone, which is approximated by a cylinder with half-height 
L = 1 — 15 kpc and radius i? = 20 kpc, and solve the equation for the 
steady-state case, where df /dt = 0. 

The first term on the right-hand side of the transport equation is the dif- 
fusion term, which accounts for the propagation through the tangled Galac- 
tic magnetic field. The diffusion coefficient K{T, r) is assumed to be constant 
throughout the diffusion zone and is parameterized by: 

K{T) = Ko(3n^ , (14) 

where P = v/c is the velocity in units of the speed of light c, and TZ is 
the rigidity of the particle, which is defined as the momentum in GeV per 
unit charge, TZ = p[GeV]/Z . The normalization Kq and the spectral index 
6 of the diffusion coefficient are related to the properties of the interstellar 
medium and can be determined from flux measurements of other cosmic- 
ray species, mainly from the Boron-to-Carbon (B/C) ratio [l5]. The sec- 
ond term accounts for energy losses due to inverse Compton scattering on 
starlight or the cosmic microwave background, synchrotron radiation and 



,2dJ_ 
dE 



9 



ionization. The third term is the convection term, which accounts for the 
drift of charged particles away from the disk induced by the Milky Way's 
Galactic wind. It has axial direction and is also assumed to be constant 
inside the diffusion region: Vc(r) = Vc sign(2;) k. The fourth term accounts 
for antimatter annihilation with rate Fann) when it interacts with ordinary 
matter in the Galactic disk, which is assumed to be an "infinitely" thin disk 
with half-height h = 100 pc. Lastly, Q{T, f) is the source term of positrons 
or antiprotons which is given by 

mnuTDM ah 

where dN/dE is the energy spectrum of antiparticlcs created via the decay 
of dark matter particles. In the transport equation, reacceleration effects 
and non-annihilating interactions of antimatter in the Galactic disk have 
been neglected. 

The solution of the transport equation at the Solar System, r = r©, 
2; = 0, can be formally expressed by the convolution 

1 /•Tmax dN(T') 

f{T) = / dT'G{T,T')—^ , (16) 

"ZDMTDM Jo dT' 

where Tmax = "i-dm for the case of the positrons and T^ax = ^^idm — 'Ti^p for 
the antiprotons. The solution is thus factorized into two parts. The first 
part, given by the Green's function G(T,T'), encodes all of the informa- 
tion about the astrophysics (such as the details of the halo profile and the 
propagation of antiparticlcs in the Galaxy). The remaining part depends 
exclusively on the nature and properties of the decaying dark matter parti- 
cles, namely the mass, the lifetime and the energy spectrum of antiparticlcs 
produced in the decay. 

We will consider the dark matter lifetime to be a free parameter that 
is constrained by requiring a qualitatively good agreement of the predicted 
positron fraction with the PAMELA results. Therefore, the only uncertain- 
ties in the computation of the antimatter fluxes stem from the determination 
of the Green's function, e.g. from the uncertainties in the propagation pa- 
rameters and the halo profile. As it turns out, the uncertainties in the precise 
shape of the halo profile are not crucial for the determination of the primary 
antimatter fluxes, since the Earth receives mostly antimatter created within 
a few kpc from the Sun, where the different halo profiles axe very similar. 
On the other hand, the uncertainties in the propagation parameters can sub- 
stantially change the predictions for the antimatter fiuxes, even by as much 
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(15) 



as two orders of magnitude for the antiproton flux. The reason for this 
large uncertainty is a correlation among the diffusion parameters and the 
size of the diffusion zone. The ranges of the astrophysical parameters that 
are consistent with the B/C ratio and that produce the maximal, median 
and minimal positron and antiproton fluxes are listed in Tab. [2] and [3] [46j. 

Positrons and antiprotons have different properties regarding their prop- 
agation, and their respective transport equations can be approximated by 
different limits of Eq. (jl3p . By exploiting the cylindrical symmetry of the 
problem, it is then possible to find semi-analytical solutions to the transport 
equation in each case. Approximate interpolating functions for the Green's 
function can be found in Ref. |10j . 
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Table 2: Astrophysical parameters compatible with the B/C ratio that yield the 
minimum (M2), median (MED) and maximal (Ml) flux of positrons. 

Positron Flux. For the case of the positrons, Galactic convection and 
annihilations in the disk can be neglected in the transport equation, which 
is then simplified to: 

V • [KiT,r)VU] + A[6(T,f)A+] + Q(r,f) = , (17) 

where the rate of energy loss, b{T,f), is assumed to be a spatially constant 
function parameterized by b{T) = ^T.^ , with Tq = 1 GeV and te = 10^^ s. 

The solution to this equation is formally given by the convolution Eq. (jl6p . 
The explicit form of the Green's function is [8] 

oo 

G,+ (T,r')= Yl i?n™(T,r')Jo(Cn^)sin(^), (18) 

n,m=l 

where Jq is the zeroth-order Bessel function of the first kind, whose succes- 
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sive zeros are denoted by Cn- On the other hand, 
teTq^ 



Bnm{T, T') 



rp2 '^nm X 



(19) 



<5-l 



with 



sm 



(20) 

where Ji is the first-order Bessel function. 

The interstehar positron flux from dark matter decay is finally given by 



$p™(r) 



4vrmDMTDM 



JTlDM 



dT'G^+{T, T' 



dT' 



(21) 



The dependence of the positron flux on the diffusion model is mostly im- 
portant at low energies, where the signal lies well below the background. 

Rather than measuring the positron flux, most experiments measure the 
positron fraction, <I>g+ / (<I>g- + since most sources of systematic errors, 
such as detector acceptance or trigger efficiency, cancel out when computing 
the ratio of particle fluxes. The background to the positron flux from dark 
matter decay will be constituted by a secondary positron flux originating 
from the collision of primary protons and other nuclei on the interstellar 
medium. 

For the background fluxes of primary and secondary electrons, as well as 
secondary positrons, we use the parameterizations obtained in Ref. |7J from 
the GALPROP numerical code for cosmic-ray propagation [2] . We leave the 
normalization of the primary electron flux as a free parameter to be fltted 
in order to match the observations of the positron fraction. 
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Ko (kpc7Myr) 


L (kpc) 


Vc (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 
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Table 3: Astrophysical parameters compatible with the B/C ratio that yield the 
minimal (MIN), median (MED) and maximal (MAX) flux of antiprotons. 
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Antiproton flux. The general transport equation, Eq. (jl3p . can be simpli- 
fied by taking into account that energy losses are negligible for antiprotons. 
Therefore, the transport equation for the antiproton density, fp{T,f,t), is 
then given by: 



= ^ = V • {K{T, f) V/p) - V • (T4(f)/p-) - 2M(z)ran„/p + Q{T, r") , (22) 



where the annihilation rate, Fann) is determined by the Galactic Hydrogen 
and Helium densities and the proton-antiproton scattering cross-section, for 
which we use the parameterization by Tan and Ng [47j . 

Analogously to the positron case, the solution to the transport equa- 
tion can be expressed as a convolution of the form Eq. (jl6l) . The analytic 
expression for the Green's function reads [18] : 



GpiT.T') = E«-P(-J;^)x (23) 



Ai{T)s\nh{Si{T)L/2) ^ V' R 



where 
and 



j\'dr' J^[CiQj X (24) 



'yc{L-z')\^^.^^^(S.{L-z' 



^^'^^^ ' \k{t) ' ^^^^ ' — ' 



Ai{T) = 2/irann(T) + T4 + A:5,(T)coth^^^ , (25) 



- ^j§y+§- (26) 
The interstellar antiproton flux is then given by 



= ^'a,(T,r)^. (27) 

where Vp{T) is the antiproton velocity. However, in order to compare the 
calculated antiproton spectrum with experimental results, one also has to 
take into account the effect of solar modulation. In the force field approxi- 
mation [491 [50] the effect of solar modulation can be included by applying 



13 



the following relation between the antiproton flux at the top of the Earth's 
atmosphere and the interstellar antiproton flux [51]: 

^^-^ox)-^ 2m,ris + ri| j'^^^(^is), (28) 

where Tjs = Ttoa + </'F; with Tjs and Ttoa being the antiproton kinetic 
energies at the heliospheric boundary and at the top of the Earth's atmo- 
sphere, respectively, and (pp being the solar modulation parameter, which 
varies between 500 MV and 1.3 GV over the eleven-year solar cycle. Since 
experiments are usually undertaken near solar minimum activity, we will 
choose (j)F = 500 MV for our numerical analysis in order to compare our 
predicted flux with the collected data. For the antiprotons, we will not ex- 
amine the background from spallation, but simply require that the fluxes 
from dark matter decay lie appropriately below the measurements so as to 
be compatible with predominantly secondary antiproton production. 

4 Results 

Below, we will present our results for the cosmic-ray signatures of decaying 
neutralinos and decaying hidden gauginos. In both cases, we will start with 
an analysis of the predictions that follow when assuming that the visible 
sector is described by an exemplary point in the coannihilation region of 
the mSUGRA parameter space. This ensures a consistent cosmology in the 
visible sector and that all free parameters of the MSSM are flxed. After that, 
we will go beyond this mSUGRA scenario and discuss how the cosmic-ray 
signatures can change in more generic cases. This will include a discussion 
about cascade decay in light of the apparent double-peak structure of the 
ATIC data. 

4.1 Decaying Neutralinos 

The potentially relevant decay modes for the case Mx < M^o, where the 

Al 

lightest neutralino can decay into the hidden gaugino, are summarized in 
Tab. [H Beside three body decays, which produce fermion/anti-fermion 
pairs, we also have to take into account the decay into Higgs- and 
bosons. Throughout the analysis we will assume that the lightest neutralino 
Xi makes up the dominant part of the dark matter, p^o ~ pdm- 

mSUGRA point. As stated above, our exemplary mSUGRA model lies 
in the coannihilation region. The deflning parameters are niQ = 150 GeV, 



14 



"^1/2 = 720 GeV, ^0 = 0, tan/? = 10 and sign/i = +1. In this model the 
Ughtest neutraUno has a mass of 301 GeV and the correct rehc density to 
be dark matter, O/i^ ~ 0.104lf] As typical for models in the coannihilation 
region, the three right-handed sleptons have masses around 304 — 307 GeV, 
which is similar to the mass of the lightest neutralino. Furthermore, this 
particular mSUGRA point, which features e.g. a spin-independent cross 
section per proton of 2.6 x 10"'*^ cm'^, can be probed with the next-to-next 
generation direct DM detection experiments like XENONIT and LUX/ZEP. 



Mx [GeV] 


e-e+X 


Branching 


Ratios for 
T-r+X 


h'^X 




r^o [1026 s] 


1 


28% 


28% 


32% 


8.8% 


2.6% 


1.8 


50 


27% 


27% 


30% 


13% 


2.4% 


1.7 


100 


24% 


24% 


28% 


21% 


2.4% 


1.5 


150 


21% 


21% 


24% 


32% 


2.6% 


1.3 


200 


30% 


30% 


36% 




3.7% 





Table 4: Branching ratios for the decay of a neutralino Xi iiito a lighter 
hidden gaugino X, for different hidden gaugino masses Mx- In the visible 
sector, masses and mixing parameters are fixed by a mSUGRA scenario in 
the coannihilation region as described in the text. The lightest neutralino 
has a mass of 301 GeV. Branching ratios of three-body decays into neutri- 
nos, Xi ~^ vvX, and quarks, Xi ~^ QQ^^ are smaller than 0.3% and 0.02%, 
respectively. The two-body decay into photons, Xi ~^ l^i is one-loop sup- 
pressed and neglected. We also indicate the lifetime of the neutralino which 
gives the best fit to the data. 

The dominant branching ratios of the neutralino decay are summarized 
in Tab. HI for different masses of the hidden gaugino Q Most interestingly, 
the fraction of decays into charged leptonic final states is never below ~ 65%. 
Beside the small masses of the right-handed sleptons, the underlying reason 
is the large /i-term, /u = 865 GeV, which suppresses the mixing between the 
bino-like lightest neutralino x?, the hidden gaugino and the higgsinos like 
~ 0{Mz/^i). 

To obtain the energy distribution of gamma rays, positrons and antipro- 
tons that are produced in the neutralino decay, we use the event generator 
PYTHIA 6.4 |55j . From these spectra, the contribution to cosmic-ray fluxes 

^The mass spectrum and relic abundance were calculated with the aid of Dark- 
SUSY 5.0.4 [sg. 

*The calculations were done with FeynArts 3.4 and FormCalc 5.4 [531 154] . 
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as measurable at Earth can be derived as described in the previous sec- 
tion. Note that the hfetime of the neutrahno is always fixed by requiring 
a qualitatively good agreement with the positron fraction as measured by 
PAMELA. 




Figure 1: Positron fraction, total electron+positron, extragalactic gamma- 
ray and antiproton flux of a decaying neutralino Xi predicted for our 
exemplary mSUGRA scenario. The used branching ratios are shown in 
Tab.m The mass of the decaying neutralino is 301 GeV, the hidden gaugino 
mass varies between 1 GeV (solid), 50 GeV (dotted), 100 GeV (dashed) and 
150 GeV (dot-dashed). We used the MED propagation model. In the lower 
left plot, the grey lines indicate the flux without background. In the lower 
right plot, we only show the flux without background. 

Our results are shown in Fig. [TJ We find that in principle the model 
can account for the observed excess in the positron fraction around 10 — 
100 GeV if the hidden gaugino is light with a mass Mx ^ 50 GeV, although 
the predicted peak seems to rise too slowly to fully match the PAMELA 
data. This slow rise is due to the two-body decay into Higgs bosons, whose 
subsequent fragmentation produces rather soft positrons. From the lower 
plots of Fig. [T] it is apparent that the model is compatible with the EGRET 
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measurements of the extragalactic gamma-ray background|j However, the 
contribution to the antiproton flux can be problematic for hidden gaugino 
masses above ^ 100 GeVH Of course, the peak in the ATIC data around 
300 — 800 GeV cannot be reproduced in this setup. 

Idealized three-body decay of a heavy neutralino. For different pa- 
rameters of the underlying MSSM model, the above plots can mainly change 
in two ways. Firstly, a larger value of the u-parameter would reduce the 
branching ratio into Higgs- and bosons As a result, the rise in the 
positron fraction would be steeper, and the contribution to the antiproton 
flux smaller. Secondly, a higher mass of the decaying neutralino would shift 
the peak to higher energies, as suggested by the ATIC data. 

In Fig. [2]we show our results for the cosmic-ray fluxes in the idealized case 
where a bino-like neutralino decays only via virtual right-handed sleptons. 
This resembles scenarios with a large term and large masses for the left- 
handed sleptons. The masses of the hidden gaugino and the neutralino Xi 
are 150 GeV and 500 GeV (sohd lines) or 300 GeV and 1850 GeV (dotted 
lines), respectively. Note that the thick lines correspond to the standard 
case where the neutralino decays democratically into all three flavors. 

As expected, the rise in the positron fraction is now steeper and can easily 
accommodate the PAMELA data. Furthermore, a very heavy neutralino 
around 1.8 TeV allows to also account for the ATIC excess. In any case we 
find a clear excess in the extragalactic gamma-ray flux at energies above 
10 GeV. 

The gamma rays come mainly from r decays and bremsstrahlung, but the 
latter is a subdominant effect as long as the three-body decay into charged 
leptons is democratic. However, the decay into taus can be suppressed in 
cases where the stau mixing angle is large, since the correspondingly larger 
left-handed component of the lighter stau weakens the coupling to the bino- 
like neutralino and the hidden gaugino. For example, if the lighter stau 
is equally left- and right-handed, the three-body decay into taus would be 

^ Gamma rays with energies below ~ 10 GeV stem from the fragmentation of the Higgs 
boson whereas gamma rays at higher energies mainly come from t decay. 

^ Note that the uncertainty in the antiproton flux at Earth from dark matter decay 
can be as large as one order of magnitude in both directions [10], due to our ignorance of 
the precise propagation parameters. 

^ A concrete lower bound on the /^-parameter is extremely model-dependent. However, 
in the concrete scenario with a bino-like lightest neutralino where we take tan/3 = 10, 
ah ~ —0.1 and assume that the right-handed sleptons have a mass around 1.02 • M^o, the 
lower bound fi > 0{2AI^o) turns out to be sufficient to suppress the branching ratio into 
Higgs-bosons below 20%. 
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Figure 2: Positron fraction, total electron+positron flux and extragalactic 
gamma-ray flux for an idealized, three-body decaying bino-like neutralino. 
We neglect effects from and bosons and assume pure democratic three- 
body decay into charged lepton pairs. The masses of the neutralino and the 
hidden gaugino are 500 GeV and 150 GeV (thick solid lines) or 1850 GeV 
and 300 GeV (thick dotted lines), respectively. The thin lines show the 
predictions when the decay into the tau-channel is neglected. The mass of 
the right-handed sleptons is assumed to be by a factor 1.1 larger than the 
neutralino mass. 

suppressed by a factor of ~ 0.5. To obtain a lower bound on the predicted 
gamma-ray signal, we also show the case where the neutralino decays into 
the first two generations only (see thin lines in Fig. ^ . The gamma-ray flux 
is much smaller in this case and comes mainly from bremsstrahlung of the 
electrons produced in the three-body decay Xi ~^ e~e~^X. Note that in any 
case we obtain a tight correlation between the contributions to the positron 
flux and the extragalactic gamma-ray flux, where the latter comes mainly 
from T decays in most cases. 

Heavy bino-like neutralinos with masses above a few hundred GeV are 
problematic for cosmology, since they are typically overproduced, even when 
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coannihilation with sleptons is taken into account. At the same time, wino- 
and higgsino-hke lightest neutralinos do not exhibit the desired leptophiUc 
decayl^ However, these problems are absent if one considers scenarios where 
the hidden gaugino is heavier than the lightest neutralino, Mx > M o. 
First, due to the mixing with the bino, the interactions of the hidden gaug- 
ino are automatically "bino-like". Second, for the small mixings that we 
consider bounds from overproduction arguments are irrelevant [29]. Note 
that the results from this paragraph can also hold in that case, provided 
one exchanges the roles of the hidden gaugino and the lightest neutralino. 
However, this requires that all sparticles, apart form the lightest neutralino, 
are heavier than the hidden gaugino. Generically this will not be the case 
and the hidden gaugino will cascade-decay through the different sparticles 
down into the lightest neutralino. We will consider this in detail in the next 
subsection. 

4.2 Decaying Hidden Gauginos 

A hidden gaugino that is heavier than the lightest neutralino, Mx > M^o, 
turns out to be more appealing from the phenomenological point of view. 
In this case, the mass of the lightest neutralino can be small and of the 
order of a few 100 GeV, and the hidden gaugino automatically possesses the 
"bino-like" interactions which are desirable for the leptophilic decay. We will 
again assume that the lightest neutralino makes up most of the dark matter, 
whereas the hidden gaugino contributes only a subdominant part px <C P^o 
to the overall matter density of the Universe. Note that in this case the 
lifetime of the hidden gaugino can be as small as tx ~ 10^^ s, the current 
age of the Universe, provided that its relic abundance is small enough. For 
definiteness, we will take px = W~^p-^o throughout this section. 

We will firstly discuss the contributions from the decaying hidden gaug- 
ino to the cosmic-ray flux as predicted for our reference mSUGRA scenario. 
Secondly we will consider the multi-peak structure of a cascade-decaying 
hidden gaugino in light of the ATIC data. 

mSUGRA point. The considerations in this paragraph are again based 
on the mSUGRA scenario described above. Depending on the mass of 
the hidden gaugino, its decay can produce fermions, neutralinos, charginos, 
Higgs- and gauge bosons as depicted in Tab. [H The corresponding branch- 

*Winos only couple to left-handed sleptons, which are typically heavier than the right- 
handed ones, whereas higgsinos can easily decay into the Higgs boson. 
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ing ratios are summarized in Tab. \5\ where we do not show the subsequent 
decays of the neutrahnos X2 3 4 charginos simphcity0 



Mx [GeV] 


I'D 


Branching Ratios for X - 
11 qq h^Xi 


X? 


w^xf 


TX [10^3 s] 


600 


1.8% 


98.2% 




0.1% 


0.0% 




1.1 


700 


5.6% 


92.9% 




0.6% 


o.c 


)% 


0.9% 




800 


5.6% 


84.6% 




3.5% 


0.^ 


!% 


6.1% 


1.0 
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0.7% 


49.8% 




17.3% 


1.^ 


!% 


31.0% 
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15.3% 


53.7% 




10.7% 
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)% 
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0.8 
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14.1% 


81.1% 




1.4% 


l.C 


)% 


2.4% 
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13.3% 


76.8% 




2.7% 


2.1 


i% 


4.7% 


0.7 


1400 


13.2% 


74.1% 


1.6% 


2.9% 
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5.4% 




1600 


12.5% 


68.5% 


8.4% 


2.7% 


2.7% 


5.2% 





Table 5: Branching ratios of the dominant decay modes of a hidden gaugino 
that is cascade-decaying into the MSSM particle zoo. The underlying sce- 
nario is our chosen mSUGRA reference point as described in the text. Neu- 
trinos and charged leptons decay essentially democratically into the three 
different flavors. We also indicate the lifetime of the hidden gaugino that 
gives the best fit to the PAMELA data, assuming an energy density of 
Px = 10"^pDM- 

As apparent from Tab.[5l the decay into charged lepton/slepton pairs is 
dominant in the whole mass range Mx — 600 — 1600 GeV that we consider. 
The decay into quarks is suppressed by the large squark masses, niq > 
1.1 TeV, whereas decay into h^, and bosons is mainly suppressed 
by the small mixing between higgsinos and the hidden gaugino. However, 
this mixing can become enhanced when the masses of the higgsinos become 
comparable to the mass of the hidden gaugino, which happens around Mx ~ 
870 GeV. 

Our results for the cosmic-ray fluxes are shown in Fig. [3] for hidden 
gaugino masses between 600 GeV and 1200 GeV, where we adjusted the 
lifetime of the hidden gaugino to fit the PAMELA data. For all considered 
masses of the hidden gaugino the predictions for the positron fraction are in 
qualitatively good agreement with the PAMELA data. At the same time, 

^These subsequent decays are taken into account in our calculations. We singled out the 
dominant decay modes in our reference mSUGRA model and used them in the PYTHIA 
code: X h\l X ^ Z\l X W^xt. X§ - xl - X°ih° and - 

xtZ°{28%), xth°{27%), X2W=^(36%). The decay of X2 and xt only produces leptons 
and is neglected. 
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Figure 3: Positron fraction, extragalactic gamma-ray flux, antiproton flux 
and total electron + positron flux from the decay of a hidden gaugino as 
predicted by our mSUGRA scenario. The branching ratios are shown in 
Tab. m The mass of the hidden gaugino varies between 600 GeV (sohd), 
800 GeV (dotted), 1000 GeV (dashed) and 1200 GeV (dot-dashed). 



the contribution to the antiproton flux hes well below the measurements 
and hence is safe in all cases. Furthermore, we obtain contributions to the 
extragalactic gamma-ray flux, which are mainly due to r decays. They are 
compatible with the EGRET measurements but could show up in future 
experiments as an excess above background. The total electron+positron 
flux is also compatible with the different measurements and we predict a 
sharp step at high energies. 



Multi-peak structures from cascade decays. As already evident in 
Fig. El the energy distribution of particles produced in cascade decays in gen- 
eral features several peaks. Their exact position carries information about 
the masses of the different intermediate particles. It is intriguing to spec- 
ulate that the apparent double-peak structure of the ATIC data originates 
from cascade-decaying particles |12j . 

In the case of the decaying hidden gaugino, the energy spectrum of 
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positrons in general possesses two pronounced peaks. These peaks stem 
from decays with intermediate selectrons. To simphfy the discussion, we will 
neglect decay modes that produce h^, and bosons or left-handed 
sleptons, and we will assume democratic decay into all three flavors. If 
we furthermore assume approximate mass degeneracy for the three right- 
handed sleptons, we are left with only three free parameters: the mass of 
the hidden gaugino Mx , the mass of the lightest neutralino M^o , and the 
mass scale of the right-handed sleptons Mr . 




100 200 300 400 500 600 700 
Energy [GeV] 



Figure 4: Energy spectrum of positrons from an idealized cascade-decaying 
hidden gaugino. Only two-body decay into right-handed slepton/lepton 
pairs is taken into account. The slepton subsequently decays into the light- 
est neutralino. The spectrum exhibits two pronounced peaks, which we 
denote by and Ei. We show plots for a lightest neutralino with mass 
150 GeV (solid) and with ITeV (dashed). The position of the peaks is fixed 
to Eh = 700 GeV and Ei = 200 GeV, as suggested by the ATIC data. The 
masses of the right-handed sleptons follow then from Eq. (|29p . We also in- 
dicate the part of the positrons that comes solely from the tau/stau decay 
channel (blue area). 

In Fig. S] we show the corresponding energy spectrum of positrons for 
two different sets of particle masses. The two pronounced peaks are denoted 
by Efi and Ei. Fixing the neutralino mass and the position of the peaks 
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determines the slepton and hidden gaugino masses according to 

Mx = E,, + ^El + Ml . 

As a simple attempt to fit the ATIC data with an ideahzed cascade- 
decaying hidden gaugino, we take the values Eh = 700 GeV and Ei = 
200 GeV. After this, a neutralino mass of M^o = 150 GeV (1000 GeV) im- 
plies a slepton mass of = 177 GeV (1117 GeV) and a hidden gaugino 
mass of Mx = 1422 GeV (2018 GeV). 
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Figure 5: Positron fraction and total electron+positron flux for an idealized 
cascade-decaying hidden gaugino. Like in Fig. [U only leptonic decay modes 
are taken into account. We assume democratic decay into the three right- 
handed charged sleptons (thick lines), or into only the selectron and smuon 
(thin lines). The mass of the lightest neutralino varies between 150 GeV 
(solid) and 1 TeV (dotted) . We show the plots for the propagation models 
MED (upper plots) and M2 (lower plots) of Ref. 06] (see Tab. [2]). 

The resulting cosmic-ray fluxes for the two neutralino masses are shown 
in Fig. [5] (upper plots, thick lines), where we used the MED propagation 
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model. Interestingly, the predicted fluxes are practically the same for the 
two cases, although the bump structure is slightly enhanced in case of the 
heavier neutralino. The lower plots are based on the M2 propagation model 
as an exemplary model with a thin diffusion zone, L = 1 kpc (as opposed to 
L = 4 kpc in case of the MED model). Since in these scenarios electrons and 
positrons are more likely to escape the diffusion zone before having lost much 
of their energy, the observable spectrum becomes steeper. This effect can 
improve the agreement with the sharp features of the ATIC data. However, 
note that propagation models with thin diffusion zones have problems with 
the correct prediction of the flux of unstable isotopes like e.g. ^^Be or 

Up to now we have assumed a vanishing stau mixing angle and demo- 
cratic decay into all three flavors. However, if the stau mixing angle is large 
the decay mode into tau/stau pairs would be suppressed, as discussed above. 
In Fig.[4]we indicated the part of the positrons that comes from the stau/tau 
channel in the case of democratic decay (blue region). A suppression of this 
channel can lead to a relative enhancement of the two peak structure of the 
cascade decay. This effect is shown by the thin lines in Fig. El where we only 
took into account the decay modes into muon/smuon and electron/selectron 
pairs. As expected, the peaks at high energy become more pronounced, and 
the spectrum becomes harder at low energies, although the effect is not 
dramatic. 

5 Conclusions 

In this work we have shown that a simple extension of the MSSM by an 
additional hidden abelian gauge group U{l)x, which kinetically mixes with 
the hypercharge U{1)y, can account for the observed PAMELA excess if 
the kinetic mixing parameter is in the range of x ~ lO"^^'^ --^''). We also 
shortly discussed possible origins of such a tiny mixing in scenarios with 
warped extra dimensions. Depending on the masses, either the visible sector 
neutralino or the hidden gaugino becomes unstable and subject to decay. We 
have demonstrated that this decay is dominated by leptonic modes in certain 
parameter regions of the MSSM where the sleptons are light (see Tab. [Hand 
Tab.E]). We found that a decaying hidden gaugino with a mass around 600 — 

^^Furthermore, the change of the propagation model in principle also changes the pre- 
dictions for the background of secondary positrons, which were calculated for the above 
Ml, M2 and MED model in Ref. [5S]. However, in Fig. [S]we used the background from 
Ref. [2] since the backgrounds obtained in Ref. [58] are meant to be extreme cases. 



24 



1200 GeV can naturally explain the observed excess in the positron fraction 
without overproducing antiprotons (see Fig. [3l). Our considerations suggest 
a preference for supersymmetric models with relatively light sleptons, as 
e.g. realized in mSUGRA models which lie in the coannihilation region. In 
any case, we predict a contribution to the extragalactic gamma-ray flux, 
which mainly stems from tau decays, and which should be observable in 
future experiments like the Fermi Gamma-ray Space Telescope. We also 
demonstrated that it is difficult to accommodate the sharp double-peak 
structure in the ATIC data within our model and with standard propagation 
models. 
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